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•  • 1 J  /  o  r 


Introduction 


It  is  often  easy  to  secure  estimates  for  the  highest  and  lowest  temperatures 
to  be  expected  in  a  structure,  but  it  is  much  more  difficult  to  obtain  detailed 
information  on  the  temperature  distribution  itself,  short  of  a  complete  solu¬ 
tion  of  the  problem.  It  has  been  found  possible,  and  quite  useful,  in  many 
cases,  to  determine  the  maximum  possible  thermal  stress  which  may  develop, 
solely  on  the  basis  of  the  known  bounds  on  the  temperature.  Upper  and  lower 
bounds  on  the  thermal  stresses  and  deformations  have  been  constructed  in  this 
way  for  beams  and  plates  [1]  and  for  composite  structures  [2,3,4].  Applica¬ 
tions  of  these  bounds  have  been  developed  in  the  estimate  of  errors  in  approxi¬ 
mate  calculations  [5]  and  in  the  analysis  of  thermal  reduction  of  beam  torsional 
rigidity  [6],  as  well  as,  of  course,  in  many  practical  design  estimates.  The 
present  paper  develops  similar  bounds  for  the  thermal  stresses  which  may  arise 
in  solid  or  hollow  spheres.  The  temperature  distributions  considered  are 
axisymmetric,  and  are  otherwise  either  arbitrary  or  radially  monotonic. 

Analysis 

For  a  hollow  sphere  of  inner  radius  a  and  outer  radius  b  with  an  axially  sym¬ 
metric  temperature  variation  T(r),  the  radial  stress  the  circumferential 

stresses  and  a  ^  and  the  displacement  u  are  [7]: 

v  ■  ((1  -  f*>  l  ^  -  »  i  <la) 

°ee  ■  *  2>  l  Tp2d0  +  $  +  2>  {  ^  -  a  -  6S,T)  (lb) 


u  "  2m(?-F)  {f*  l  Tp2dp  +  pr  ?  Tp2dp  +  2P\^2V)  {  Tp2dp> 

fl  r  E 

where  8  ■  — ,  p  ■  — ,  K  •  •; — ;  E  is  Young's  modulus,  u  the  shear  modulus, 

D  D  l“v 

Poisson's  ratio  and  a  the  coefficient  of  thermal  expansion. 

Assume  now  that  bounds  on  the  temperature  distribution  are  known,  i.e. 


Tm  1  T(r)  1 

Ue  can  then  write,  from  equation  (la) , 


"rr^  :  -PF  <«  -  -  ‘X^V 


2K 


and  therefore  a  bound  for  a 


rr 


p~'  3  /XM  vp3  i,v  3 
at  any  given  value  of  p  is: 


Iff  (p)  |  <  — -r  (p3~83)  -  T  ) 

1  rr'p*1  -  T-fF  3P3  llM  m' 


The  expression  modifying  (Tu  -  T  )  is  positive  and  assumes  a  maximum  at 

n  m 

p2  *  8,  and  so  in  the  interval  6  1  P  <  1J 

.  2K  1— 2 

Vl 

For  the  case  of  a  thin  shell,  x  ■  is  small,  and  equation  (4)  can  be 
written  as: 


rr1  - 


k(tm  - 


V(7  +  T  + 


7 

961 


Similarly,  we  obtain  from  equations  (lb)  and  (lc)  the  bounds: 


<  K(Tm  -  T  ) 
“Mm 


(2) 


(3) 


(4) 


(5) 


(6a) 


aaT  <  u  <  obT. 
m  -  1 


(6b) 


Monotonic  Temperature  Distributions 


When  the  temperature  increases  monotonically  with  r,  then  T  -  T(B)  and 

m 

T(l)  -  T  ,  and  the  lower  bound  on  a  can  be  shown  to  be  zero  since,  from 
M  rr 

equation  (la),  we  have: 


(T(p)-^<1  -  |J)  -  T(p)  (~j  -  1)^4 


0 


so  that 


and  as  before 


|o00i  <  KfT(l)  -  T(B)I 


(7) 


(8) 


aaT(B)  <  u  <  obT(l) 


(9) 


Stresses  at  the  Edges 

For  an  arbitrary  temperature  distribution,  the  circumferential  stresses  at  the 
inner  and  outer  edges  of  the  sphere  are,  respectively: 

o00(8)  -  t3  l  Tpzdp  -  (1  -  B3)T(8» 


aee(1) 


Tp2dp  -  (1  -  63)T(1)> 


Hence, 


a00(i)  -  o00(B)  +  K[T(B)  -  T(l)] 


(10) 


Using  the  relation  in  equation  (8),  we  conclude  that,  for  a  monotonically 
increasing  temperature  distribution,  Oqq  (1)  is  a  compressive  stress  and 

(6)  is  a  tensile  stress.  The  reverse  holds  for  a  monotonically  decreasing 
temperature  distribution.  Note  that  T(6)  and  T(l)  cannot  be  equal  unless  T 
is  constant  throughout. 

Solid  Sphere 

For  a  solid  sphere  a  ■  0,  and  equations  (1)  reduce  to: 

°rr  "  2K{  l  Tp2dp  '  JT  l  Tp2dp} 


0ee  “  aH  "  K  {  p*  l  Tp2dD  +  2  l  Tp2dp  ■  T}  (11) 


u  "  2P  l  Tp2dp  +  f0  Tp2dp} 

Again,  with  T  <  T(r)  <  T„,  we  can  write  the  bounds  as: 
m  —  —  n 

\oJ  <  (2/3)K[Tm  -  TJ 


I  "09 1  i  KtTM  -  T„1 


oT 


m 


<  u  <  oT„ 
—  —  m 


and  for  a  monotonically  increasing  temperature  distribution: 


0  -  °rr  -  "  T(0) ]  <  |k[T(1)  - 


-KtT(l)  -  T(0)]  <  a0Q  <  |  K[T(1)  -  T(0)] 


aT(O)  <  u  <  aT(l) 


It  can  be  noted  that  for  the  solid  sphere 

°rr(0)  “  °qq<0)  "  2K{|  Tp2dp  - 

and  the  equilibrium  condition  [7] 

d8  .  do 

TT  ‘  r  (°rr  '  °e8)  ?ielda  "dT  '  0 

when  r  *  0.  However,  for  the  hollow  sphere,  <*rr(6)  5  0  and 

1 

lim  oflfl(S)  -  K{  3  /  Tp2dp  -  T(0)> 

0*0  06  ° 
a 

y 

At 


T(0)  ] 


consequently. 


becomes  infinite  as  0  0. 
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